Quantum phases of Fermi-Fermi mixtures in optical lattices 
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The ground state phase diagram of Fermi- Fermi mixtures in optical lattices is analyzed as a func- 
tion of interaction strength, population imbalance, filling fraction and tunnehng parameters. It is 
shown that population imbalanced Fermi-Fermi mixtures reduce to strongly interacting Bose-Fermi 
mixtures in the molecular limit, in sharp contrast to homogeneous or harmonically trapped systems 
where the resulting Bose-Fermi mixture is weakly interacting. Furthermore, insulating phases are 
found in optical lattices of Fermi-Fermi mixtures in addition to the standard phase-separated or 
coexisting superfluid/excess fermion phases found in homogeneous systems. The insulating states 
can be a molecular Bose-Mott insulator (BMI), a Fermi- Pauli insulator (FPI), a phase-separated 
BMI/FPI mixture or a Bose-Fermi checkerboard (BFC). 

PACS numbers: 03.75.Hh, 03.75.Ss, 05.30.Fk 



I. INTRODUCTION 

In recent experiments, superfluidity of ultracold ^Li 
Fermi atoms with population imbalance was investigated 
in gaussian traps from the Bardeen-Cooper-Schrieffer 
(BCS) to the Bose-Einstein condensation (BEC) limit [l], 
SjSIJI- In contrast to the crossoverphysics found in pop- 
ulation balanced systems @j S S HI j a phase transition 
from superfluid to normal phase as well as phase sepa- 
ration were observed [lO, [lH • More recently, experi- 
mental evidence for superfluid and insulating phases has 
also been reported with ultracold ^Li atoms in optical 
lattices without population imbalance !l2j. after over- 
coming some earlier difficulties of producing Fermi su- 
perfluids from an atomic Fermi gas or from molecules of 
Fermi atoms in optical lattices [H, [13, [H, [l^ [l^l ■ 

Before and after experiments, the possibility of in- 
teresting phases in population imbalanced ultracold 
fermions have attracted intense theoretical interest [9l.[l0l. 
[n!, liiM [11 il, [IIM ilM ii . These works 
were focused on homogeneous or harmonically trapped 
Fermi superfluids, where the number of external control 
parameters is limited. For instance, in addition to popu- 
lation imbalance P and scattering parameter a^., several 
other experimental parameters can be controlled in opti- 
cal lattices such as the tunneling matrix element t^- be- 
tween adjacent lattice sites, the on-site atom-atom inter- 
actions g, the flUing fraction F, the lattice dimensionality 
T> and the tunneling anisotropy rj — ^x/^T UM- Therefore, 
optical lattices may permit a systematic investigation of 
phase diagrams and correlation effects in Fermi systems 
as a function of P, F, V and rj, which has not been 
possible in any other atomic, nuclear and/or condensed 
matter systems. 

Arguably, one-species or two-species fermions loaded 
into optical lattices are one of the next frontiers in ul- 
tracold atom research because of its greater tunability. 
The problem of one-species population balanced fermions 
in lattices has already been studied in the BCS and 



BEC limits [6|, |30|, |3l|, while recent works on popula- 
tion and/or mass imbalanced Fermi- Fermi mixtures in 
lattices were limited only to the BCS regime [3^[3^[33|. 
In this manuscript, we extend our previous work describ- 
ing balanced fermions and their s-wave and p-wave su- 
perfluid phases in optical lattices [3l|, and extend our 
recent Letter [s^ describing superfluid and insulating 
phases of Fermi-Fermi mixtures in optical lattices. We 
discuss specifically the population balanced and imbal- 
anced Fermi-Fermi mixtures of one-species (e.g. ^Li or 
^•^K only) and of two-species (e.g. ^Li and '^''K; ^Li and 
^''Sr; or '*°K and *''Sr), as a function of g, P and F for 
fixed values of rj in optical lattices and throughout the 
BCS to BEC evolution. We also pay special attention 
to the emergence of insulating phases in the strong at- 
traction (molecular) limit. The existence of insulating 
phases in optical lattices should be contrasted with its 
absence in homo gen eous or harmonically trapped sys- 
tems [m, [23I l28j. Furthermore, the present exten- 
tion to two-species mixtures is timely due to very re- 
cent experimental work on ^Li and "^^K mixtures (Sa . [37j 
in harmonic traps, which have opened the possibility of 
trapping these mixtures in optical lattices as well. 

Our main results are as follows. Using an attractive 
Fermi-Hubbard Hamiltonian to describe one-species or 
two-species mixtures, we obtain the ground state phase 
diagram containing normal, phase-separated and coex- 
isting superfluid/excess fermions, and insulating regions. 
We show that population imbalanced Fermi-Fermi mix- 
tures reduce to strongly interacting (repulsive) Bose- 
Fermi mixtures in the molecular limit, in sharp con- 
trast with homogenous systems where the resul ting Bose- 
Fermi mixtures are weakly interacting [l9l . [23l . l28j . This 
result is a direct manifestation of the Pauli exclusion 
principle in the lattice case, since each Bose molecule 
consists of two fermions, and more than one identical 
fermion on the same lattice site is not allowed for sin- 
gle band systems. This effect together with the Hartree 
energy shift lead to a filling dependent condensate frac- 
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tion and to sound velocities which do not approach to 
zero in the strong attraction hmit, in contrast with the 
homogeneous case [l!]. Lastly, several insulating phases 
appear in the molecular limit depending on the filling 
fraction and population imbalance. For instance, we find 
a molecular Bose-Mott insulator (superfluid) when the 
molecular filling fraction is equal to (less than) one for 
a population balanced gas where the fermion filling frac- 
tions are identical. When the filling fraction of one type 
of fermion is one and the filling fraction of the other is 
one-half (corresponding to molecular boson and excess 
fermion filling fractions of one- half), we also find either 
a phase-separated state consisting of a Fermi-Pauli in- 
sulator (FPI) of the excess fermions and a molecular 
Bose-Mott insulator (BMI) or a Bose- Fermi checkerboard 
(BFC) phase depending on the tunneling anisotropy 77. 

The rest of the manuscript is organized as follows. Af- 
ter introducing the Hamiltonian in SeclTTJ first we derive 
the saddle point self-consistency equations, and then ana- 
lyze the saddle point phase diagrams at zero temperature 
in Sec. mil In Sec. lIV| we discuss Gaussian fluctuations at 
zero temperature to obtain the low energy collective ex- 
citations, and near the critical temperature to derive the 
time-dependent Ginzburg-Landau (TDGL) equations. In 
Sec. |Vl we derive an effective Bose-Fermi action in the 
strong attraction limit, and describe the emergence of in- 
sulating BMI, BFC and separated BMI/FPI phases, and 
in Sec. IVII we suggest some experiments to detect these 
insulating phases. A brief summary of our conclusions is 
given in Sec. IVIII Lastly, in App. |^ and [Bl we present 
the elements of the inverse fluctuation matrix, and their 
low frequency and long wavelength expansion coefficients 
at zero and finite temperatures, respectively. 



has the form Va{r) = ya,aJ2j={x,y.z}^^'^^i'^'^j)^ where 
Vo.cr is proportional to the laser intensity, k, ~ lnjX is 
the wavelength of the laser such that — A/2 is the 
size of the lattice spacing, and Vj corresponds to the jth 
component of r. This potential describes a cubic op- 
tical lattice since the amplitude and the period of the 
potential is the same in all three orthogonal directions. 
Furthermore, we assume short-range interactions and set 
C/(r, r') = [/(5(r — r'), where C/ > is the strength of the 
attractive interactions. 

For this Hamiltonian, the single particle eigenfunc- 
tions are Bloch wave functions, leading to a set of Wan- 
nier functions that are localized on the individual lattice 
sites [Hj. Therefore, we expand the creation and anni- 
hilation field operators in the basis set of these Wannier 
functions Waij: — Ti) of the lowest energy states of the 
optical potential near their minima, such that '(/'^(r) = 
Ei W<,(r-ri)a,,„ and ^/'^(r) = Yu^ W^^(i'-ri)4<„ where 

alo- (oi.CT) site operator creates (annihilates) a fermion at 
lattice site i with a pseudo-spin a. Here, Ti is position 
of the lattice site i. Multi-bands are important when the 
interaction energy involved in the system is comparable 
to the excitation energies to the higher bands, and these 
effects can be easily incorporated into our theory. Notice 
that, for a complete set of Wannier functions such that 
Wa(j^ - ri)iy*(r' - ri) = (5(r - r'), the field operators 

\i\ja(v),i\}\,{x:')\ ~ — v' ) as well as site operators 

[flj o-, g./] — Sa.a'Si.j obey the Fermi anti-commutation 
rules, where S(r) is the delta function and 6ij is the 
Kronecker-delta. 

This expansion reduces the continuum Hamiltonian 
given in Eq. ([T|) to a site Hamiltonian 



II. LATTICE HAMILTONIAN 

To describe mixtures of fermions loaded into optical 
lattices, we start with the following continuum Hamilto- 
nian in real space 



2™, 



V'.(r) 



dr I (ir'n|(r)C/(r, r')n|(r'). 



(1) 



where V'i(r) [^'^(r)] held operator creates (annihilates) 
a fermion at position r with a pseudo-spin a, mass rria 
and chemical potential ficr- Here, Va{v) is the optical 
lattice potential, n„{v) ~ ■!/'J-(i")V'ct(i') is the density op- 
erator, and U{r, r') describes the density-density inter- 
action between fermions. We allow fermions to be of 
different species through TOq-, and to have different popu- 
lations controlled by independent /lo-, where the pseudo- 
spin cr labels the trapped hyperfine states of a given 
species of fermions, or labels different types of fermions 
in a two-species mixture. The optical lattice potential 



i,j,k,l 



(T ' ''i.i:<T 



(2) 



JdrW;ir ~ ri)[-VV(2r 



where i^-' = 
Vcr(r)]VFCT(r ~ Tj) is the tunneling amplitude be- 
tween sites i and j, Uij-k,i = U J drW^{r — ri)T4^|(r — 
rj)T4^*(r — rk)Vl^x('' ~ ^i) is the elements of strength 
of attractive density-density interactions between 
sites and {fc, Z}, and nij;cr = ^la^j.cr- this 

manuscript, we allow tunneling and interactions up 
to the nearest-neighbors, and choose ij;-' = tcrSi_j±i 
and Ui^j-u = gStjdj^kSk,i + hSij5j^k±iSk,i, respectively. 
Here, g > and h > are strengths of the on-site 
and nearest-neighbor interactions, respectively. This 
reduces the general site Hamiltonian given in Eq. ^ to 
a nearest-neighbor site Hamiltonian 



H^ = - E ^ 



E 



5'E"*'^;T"*,i;i - ^lE^'-^^T^ija' (3) 
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where («, j) restricts sums to nearest-neighbors only such 
that j — i ± 1. This is the nearest- neighbor Fermi- 
Hubbard Hamiltonian for a simple cubic lattice. 

Finally, we Fourier transform the site op- 
erators to the momentum space ones, such 
that a^,^ = il/VM)J2icO-k,<re~'^''' and a|^^ = 
(1/Vm) X^k '^k cr^*'^ where M is the number of 
lattice sites, and ^ (ak.o-) operator creates (annihi- 
lates) a fermion at momentum k with a pseudo-spin a. 
This leads to the general momentum space Hamiltonian 



with an on-site attractive interaction .9 > 0. Here, aj^ ^ is 
the fermion creation and ^ = '^k+q/2,T^-k+q/2,j, 

r:(k) 

is the pair creation operator. In addition, ^k,CT = ek,cr — /^o- 
describes the nearest-neighbor tight-binding dispersion, 
Ek.cr = 2icr6'k with 



j = {x.,y,z} 



[1 - cos(fcjac)]. 



(6) 



Hm — Ck,cr 



'^k.k;, 



^ nki,k2;Tt^(ki - k2)nk3,k4;i, (4) 



ki+k3=k2+k4 



where ^k.o- = ek,<T — fJ-a describes the nearest- 
neighbor tight-binding dispersion ek.o- — 



j={x,y,z} 



cos{kjac) where kj is the jth 



where Jla- — Ha ~ Vn.a and Vn.a is a possible Hartree en- 
ergy shift. Notice that, we allow fermions to be of differ- 
ent species through 1^, and to have different populations 
controlled by independent pto-- 

Furthermore, the momentum space sums in Eq. ([5]) 
and throughout this manuscript are evaluated as follows. 
For large and translationally invariant systems consid- 
ered here, we can use a continuum approximation to sum 
over the discrete k levels, and for a three dimensional lat- 
tice write 



component of k, and 



■LLh ui jv, aim tic is thc lattice spac- 

Here, nk,k';<T = flk.^ak'.a, and we used ^F{k.^,ky,k^) = f dkxdky 
k.(ri-rj) ^ 2V. / .cos^k.a^. There- k J-T^/a^ (W^ 



mg. 



■ik- (r;— Tj ) 



{x,y,z} 



cos{kjac 



dkrf dJzy 



F{kx, ky,kz), 



fore, in order to achieve thc momentum space 
Hamiltonian given in Eq. (|4|) from the real space 
Hamiltonian given in Eq. Ilj, we performed 
^,(r) = (l/VM)E,,k«k,.VK,(r - ri)e-'''-' and 

TAjCr) = (l/VM)E,,kaLw/;(r - ri)e^-'^^, which 
corresponds to the total transformation from the real 
space field operators to the momentum space ones. 

The second term in Eq. (g]) U{k - k') = 

is the Fourier transform 
of the interaction term and is given by [/(k — k') = 
g + '2hJ2j={x y 2} cos[(fcj — k'^)ac\. In momentum space, 
the simplification comes from the separability of the in- 
teraction potential, such that J7(k— k') — (7rs(k)rs(k')-f 
/i^^-^^g rj(k)rj(k') where the sum over c corre- 
sponds to different pairing symmetries. For a three- 
dimensional lattice, we obtain the following components: 
rs(k) — 1 for the s-wave symmetry; re(k) = [cos(fca;ac) + 
cos^kyGc) + cos{kzac)]/ for the extended s-wave sym- 
metry; r<i4(k) = [cosik^ttc) + cos{kyac) - cos{kzac)]/ V2, 
rd,2(k) = [cos(fca;ac) — cos(fcyac) + cos(fc^ac)]/%/2 and 
Trf.aCk) = [cos^k^ac) — cos{kyac) — cos{kzac)]/V^ for the 
d-wave symmetry; and rp^i(k) — \/2sm{kxac), rp^2(k) = 
-^2 sin(fcj^ac) and rp_3(k) = ^/2sm{kzac) for the p- 
wave symmetry. Notice that, symmetrized and anti- 
symmetrized combinations of cosines and sines help to 
exploit the symmetry of the lattice. In this manuscript, 
we set nearest-neighbor attraction to zero (h = 0), and 
consider only the s-wave on-site attractions {g ^ 0). 

Thus, to describe mixtures of fermions loaded into op- 
tical lattices, we use the s-wave single-band Hamiltonian 



= M I F{kx,ky,kz){7) 



H 



E^k 



9 E 

k,k',q 



fell -^k',q, 



(5) 



where F(kx,ky,kz) is a generic function of k^, ky and 
kz, L = Mac is the size of the lattice, and ki = kiUc 
is dimensionless. This continuum approximation is valid 
for large systems only where M ^ 1. Notice also that, 
flc provides a natural cutoff to the k space integrations 
in the lattice case. 

Unlike recent works on fermion pairing in optical lat- 
tices which were restricted to the BCS limit ^32, H^, we 
discuss next the evolution of superffuidity from the BCS 
to the BEC regime 3l| and the emergence of insulating 
phases [ssj. We ignore multi-band effects since a single- 
band Hamiltonian may be sufficient to describe the evo- 
lution from BCS to BEC physics in optical lattices [39| . 
Multi-bands are important when the fermion filling frac- 
tion is higher than one and/or the optical lattice is not in 
the tight-binding regime, but these effects can be easily 
incorporated into our theory. 



III. SADDLE POINT APPROXIMATION 

In this manuscript, we use the functional integral for- 
malism described in Ref. 0, B IH, HI]. The general 
method follows the same prescription as for the homo- 
geneous case dl], and we do not repeat the same anal- 
ysis for the lattice case discussed here. The main differ- 
ences between the momentum space expressions for the 
homogeneous and the lattice case come from the peri- 
odic dispersion relation for the lattice system, as given 
in Eq. ®. 
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A. Saddle Point Self-Consistency Equations 

For the Hamiltonian given in Eq. ([5]), the saddle point 
order parameter equation parallels that of the homoge- 
neous system 28], and leads to 



(a) 



g A/4- 



1 - /(i?k,l) - /(i?k.2) 

2-Ek + 



|r.(k)|^ 



(8) 



where M is the number of lattice sites, f{x) 
l/[exp(x/T) + 1] is the Fermi function. 



lA, 



|2U/2 



- 7s ^k - 



(9) 



is the quasiparticle energy when 71 = 1 or the nega- 
tive of the quasihole energy when 72 — —1, and -Ek,± = 
(i?k,i ± -Bk,2)/2. Here, Ak = Aors(k) is the order pa- 
rameter and ^k,± = ek,± — M±7 where ek,± — '2t±9]^ with 
t± = {t-[ ± ti)/2 and Ji± = (pt| ± Jj-i)/'2.. Notice that the 
symmetry between quasiparticles and quasiholes is bro- 
ken when ^k,- 7^ 0. The order parameter equation has to 
be solved self-consistently and leads to number equations 

iVt = E [l"kl V(i?k,i) + K'kl'/(-i?k,2)] , (10) 

k 

=Y. [I^kl V(i?k,2) + |^^kP/(-i?k,i)] , (11) 
k 

which are derived the same way as in the homogeneous 
case pj. Here, |ukP = (1 fk.+/-Ek,+)/2 and |wkP = 
(1 — ^k.+/£'k,+)/2. The number of a-type fermions per 
lattice site is given by 



0<n. = ^<l. 



(12) 



Thus, when ^ n^, we need to solve all three self- 
consistency equations, since population imbalance is 
achieved when either i?k.i or £'k,2 is negative in some 
regions of momentum space, as discussed next. 



B. Saddle Point Phase Diagrams 
at Zero Temperature 

To obtain ground state phase diagrams, we solve 
Eqs. ([5]), pj}]) and ([TT|) as a function of interaction 
strength population imbalance and total filling frac- 
tion 



„ nt — Til 



< F 



< 1, 



(13) 
(14) 



respectively [35|, and consider two sets of tunneling ra- 
tios rj — ti/t-[. The case of 77 = 1 is shown in Fig. [1] and 
the case of 77 = 0.15 is shown in Fig. [51 While rj — 1 cor- 
responds to a one-species (two-hyperfine-state) mixture 




FIG. 1: Phase diagrams for a one-species (®Li or ^"K) mix- 
ture of two-hyperfine states with 77 = 1: (a) versus n^, 
and (b) P versus F, for g = 5t+ and g = 10t+. The normal 
regions (outside the 'football' boundaries) and coexistence of 
superfluidity with excess fermions (CSE) and/or phase sep- 
aration (PS) (inside the 'football' boundaries) are indicated. 
The CSE/PS (normal) region expands (shrinks) with increas- 
ing fermion attraction. 



such as ^Li or ""^K, 77 = 0.15 corresponds to a two-species 
mixture (one-hyperfine-state of each type of atom) such 
as <^Li and 4°K; ^Li and ^^Sr; or and *^Sr. 

Generally, lines AB (0 < 7i| < 1; = 0) and ED 
(7i| = 0; < 7i| < 1) in Figs. [1] and [U correspond 
to normal a-type Fermi metal for all interactions, while 
points B (tt-i — l,ni — 0) and D (n-^ = 0, = 1) cor- 
respond to a Fermi-Pauli (band) insulator since there is 
only one type of fermion in a fully occupied band. Thus, 
the only option for additional fermions (t in case B and 
J, in case D) is to fill higher energy bands if the optical 
potential supports it, otherwise the extra fermions are 
not trapped. For the case where no additional bands are 
occupied, we label the corresponding phase diagram re- 
gions as 'Inaccessible' in Figs.[TJb) and^b), since either 
?T,| > 1 or 77,| > 1 in these regions. This phase bound- 
ary between the 'Normal' and the 'Inaccessible' phase is 
given by P = 1/F - 1 for P > 0, and P = 1 - 1/F for 
P < 0. 

The population balanced line AC ends at the special 
point C, where n'\ = ni — 1. This point is a Fermi-Pauli 
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(band) insulator for weak attraction since both fermion 
bands are fully occupied. Furthermore, for very weak 
attraction, lines BC (n| — 1,0 < < 1) and CD (0 < 
n| < 1, n| = 1) correspond essentially to a fully polarized 
ferromagnetic metal (or half-metal), where only the type 
of fermions with filling fraction less than one can move 
around. 

In the phase diagrams shown in Figs. [1] and [2l we indi- 
cate the regions of normal (N) phase where |Ao| = 0, and 
group together the regions of coexistence of superfluidity 
and excess fermions (CSE) and/or phase separation (PS), 
where |Ao| ^ 0. When F ^ 1, the phase diagrams are 
similar to the homogenous case [Ifl, Hll, 123, i28i| , and the P 
versus F phase diagram is symmetric for equal tunnelings 
as shown in Fig. \V[h) , and is asymmetric for unequal tun- 
nelings having a smaller normal region when the lighter 
band mass fermions are in excess as shown in Fig. [2Ib). 
Here, we do not discuss separately the CSE and PS re- 
gions since they have already been discussed in homoge- 
neous and harmonically trapped systems [13, Ell, [H, [ll] , 
and experimentally observed 0,I3,MI3|j but we make two 
remarks. 




FIG. 2: Phase diagrams for a two-species (^Li and ''"K) mix- 
ture of two-hyperfine states with rj = 0.15: (a) n| versus n|, 
and (b) P versus F, for g — 5t+ and g — 10t+. The normal 
regions (outside the 'football' boundaries) and coexistence of 
superfluidity with excess fermions (CSE) and/or phase sep- 
aration (PS) (inside the 'football' boundaries) are indicated. 
The CSE/PS (normal) region expands (shrinks) with increas- 
ing fermion attraction. 



First, the lattice phase diagram is also different from 
the homogenous systems' in relation to the to polo gical 
quantum phase transitions discussed in Ref. [23|, l28l |. 
These phases are characterized by the number of zero- 
energy surfaces of E'k,! and £'k,2 such that (I) -E'k,i(2) tias 
no zeros and -Ek,2(i) has only one, (II) i?k,i(2) has no zeros 
and £'k.2(i) has two zeros, and (III) i?k,i(2) and -Ek,2(i) 
has no zeros and are always positive corresponding to the 
P — limit. The topological phases characterized by the 
number (I and II) of simply connected zero-energy sur- 
faces of i?k,s may lie in the stable region of CSE, unlike in 
the homogeneous case where the topological phase II al- 
ways lies in the phase separated region for all parameter 
space [23l[28|. 

To understand the topology of the uniform superfluid 
phase, we define B± = p(|/(4tf) ± /i|/(4i|). For the s- 
wave symmetry considered, the zeros of i?k,s occur at 
real momenta 6*^ = B+±[B'i-\Ao\'^ / {At^ti)]'^/'^ provided 
that |Aop < 4\B_\^tfti for B+ > 0, and |Aop < -fi^V-i 
for B+ < 0. The transition from phase II to I occurs 
when 0^ 0, indicating a change in topology in the 
lowest quasiparticle band. The topology is of phase I if 
B+ < 0, and phase II is possible only when _B+ > 0, and 
therefore phase II (I) always appears in the BCS (BEC) 
side of the phase diagram. In the particular case when 
i| — ti, the topology is of phase I if ^+ < 0, and phase 
II is possible only when Jl^ > 0. 

Notice that, pt+ goes below the bottom of particle 
(hole) band as g increases, since the Cooper pairs are 
formed from particles (holes) when F < 0.5 {F > 
0.5) dJ. This leads to a large and negative (positive) 
/i+ for particle (hole) pairs in the BEC limit, thus lead- 
ing to a uniform superfluid with topological phase I (II) . 
Therefore, the topology of the entire superfluid phase is 
expected to be of phase II for F > 0.5. However, when 
i| — ti, the boundary between phase II and phase I lies 
around g « 13t+ for F = 0.1 and 5 - 30t+ for F = 0.3. 

Similar topological phase transitions in nonzero an- 
gular momentum superfluids have been discussed pre- 
viously in the literature in connection with ^He p- 
wave phases 0, and cuprate d-wave superconduc- 
tors [4l|, and more recently a connection to ultracold 
fermions which exhibit p-wave Feshbach resonances was 
made 42, 4^, j45j. However, we would like to em- 
phasize that the topological transition discussed here is 
unique, since it involves a stable s-wave superfluid, and 
may be observed for the first time in future experiments. 

Second, the phase diagram characterized by normal, 
superfluid (CSE or PS), and insulating regions may be 
explored experimentally by tuning the ratio .9/^+, total 
filling fraction i^, and population imbalance P as done in 
harmonic traps [l|, [3, S, 3 • We would like to emphasize 
that our saddle point results provide a non-perturbative 
semi-quantitative description of the normal/superfluid 
phase boundary for all couplings and tells us that the sys- 
tem is either superfluid (|Ao| 7^ 0) or normal (|Ao| = 0), 
but fails to describe the insulating phases. Thus, flrst, 
we analyze Gaussian fluctuations at zero and finite tem- 
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peratures, and then show that insulating phases emerge 
from fluctuation effects beyond the saddle point approx- 
imation. 



IV. GAUSSIAN FLUCTUATIONS 

In this section, we follow closely the formalism devel- 
oped to deal with fluctuations for homogeneous superflu- 
ids 0, S, SIH- We determine the collective modes for 
Fermi superfluids in optical lattices at zero temperature 
and the efffective equation of motion near the critical 
temperature. 



A. Gaussian Fluctuations 
at Zero Temperature 

Next, we analyze the zero temperature Gaussian fluc- 
tuations in the BCS and BEG limits for P — Q, from 
which we extract the low frequency and long wavelength 
collective excitations. For the s-wave symmetry consid- 
ered, the collective excitation spectrum is determined 
from the poles of the fluctuation matrix in much the 
same way as in homogenous systems [11]. Thus, following 
Ref. [2§|, we obtain the condition 



det 



A + C\(i\^-Dw^ 
—iBw 



iBw 
Q|q|2-i?u;2 



= 0, 



(15) 



where the expansion coefficients A, B, C, D, Q and R 
are given in App. |^ There are amplitude and phase 
branches for the collective excitations, but we focus only 
on the lowest energy Goldstone phase mode with disper- 
sion w(q) = w|q|, where 



/ AQ 



\AR + B^ 



1/2 



(16) 



is the speed of sound. First, we discuss analytically the 
sound velocity in the weak and strong attraction limits, 
and then calculate numerically the evolution in between 
these limits. 

In the weak attraction (BGS) limit when < /i+ < 
4I?<+, the coefficient that couples phase and ampli- 
tude fields vanish {B = 0), and the phase and ampli- 
tude modes are decoupled. We also obtain A = A/'(0), 
C = Q/3 = A/'(0)w|/(12I?|Aop), and D ^ R/3 ^ 
j\A(0)/(12|Aop), where V is the number of dimensions. 
Here, Af{x) — (5(^14,+ — x) is the density of states 
and Vp = A/'r(0)/A/'(0) is the Fermi velocity, where 5{x) 
is the delta function and Nr{x) = X^k iK^^^kj)^ — 
(Vi^k,j,)^]^(^k,-i- — x) is an effective 'kinetic' density of 
In addition, we find |Ao| = 2[{2Vt+f - {J1+ - 
l/(gA/'(0))] for the order parameter and 
/i+ = F/J\f{0) for the chemical potential. 

The calculation given above leads to = Vp /T>, which 
reduces to the Anderson-Bogoliubov relation when t-\ — 



states. 

2I?t+)2]i/2exp 



ti- In this limit, the sound velocity can also be written 
as = vf^ivf,i/V where Vp^^ = A/'r,o-(0)/A/'<j(0). Here, 
Na{x) = X]k^('?k(T — x) is the density of states, and 

^rA^) = Ek.(vi^k..)'<5(a,. - X). 

In the strong attraction (BEG) limit when /i+ < for 
F < 0.5 and pt+ > 4I?t+ for F > 0.5, the coefficient 
B ^ indicates that the amplitude and phase fields are 
coupled. We obtain A = 4i^(l - F)/g, B = {1- 2F)/g'^, 
D = {1- 2Ff/g^, Q = 2alt^ti/g^ and R = l/g^, where 
we used |Ao| = g[F{\ — F)Y/'^ for the order parameter 
and pt+ ^ -{g/2 + 4Pt|t|/.g)(l - 2F) - 2Vt+ for the 
average chemical potential. Notice that, /i+ — eb/2 — 
21)1+ — A'Dt-[ti{l — 2F)/g, where = —g is the binding 
energy defined by 



E 



|r.(k)p ^ 

2ek,+ — 



(17) 



and Vh — gF is the Hartree energy. Thus, we obtain 
= 8a^t|t|i^(l — F) = v-^vi for the sound velocity, 
where v„ = 2ac t„y/2F{l - F). Notice that, -y has a 
dependence in the dilute limit which is consistent with 
the homogenous result where the sound velocity depends 
on the square-root of the density. However, v saturates to 
a finite value as g increases which is sharp in contrast with 
the vanishing sound velocity of homogenous systems [1, 

Hi. 

When F < 0.5 {F > 0.5), making the identification 
that the number of particle (hole) pairs per lattice site 
is ub = F {riB = 1 — F), we obtain Ubb = 2g{l — F)al 
{Ubb — '^gF'^l) as the repulsive particle (hole) boson- 
boson interaction such that the Bogoliubov relation = 
UBBnB / (iTiBac) is recovered. Here, we also identified 
niB = (7/(4a^t|i|) as the mass of the bound pairs to be 
discussed in Sec. IIVBI Therefore, in both low and high 
filling fraction limits, we find that Ubb ~ '^9<^c- 
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FIG. 3: Sound velocity v (in units of t+ac) versus F in 
the weak attraction (WA) limit for g = 2t+ in two (hollow 
squares) and (solid squares) in three dimensions, and in the 
strong attraction (SA) limit for g = 20t+ (dotted line) in both 
two and three dimensions. Here, t = t-f = ti and P — 0. 

In Fig. [21 we show the sound velocity v for two and 
three dimensional lattices as a function of F when t = 
t| — and P = 0. In the BGS limit shown for g = 2t+, 
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V is very different for two and three dimensional lattices 
due to van Hove singularities present in their density of 
states [sO, El]. However, in the EEC hmit shown for g = 
20i+, V saturates to a finite value, which is identical for 
both two and three dimensional lattices and reproduces 
the analytical results discussed above. 

Before concluding this section, we discuss the fraction 
of condensed pairs to understand further the saturation 
of sound velocity in the BEC limit. In the absence of a 
population imbalance (P = 0), the number of condensed 
pairs per lattice site is given by [4^ 



1 



■ tanh 



(18) 



where 0(k) = Ak/(2i?k,+) is the ground state pair wave 
function. At zero temperature, while the number of con- 
densed pairs ric = 7rA/'(0)| Ao|/4 is very small in compar- 
ison to the total number of fermions in the BCS limit, 
it increases as a function of interaction strength and 
saturates in the BEC limit to ric = ^(1 — -F), leading 
to the number of condensed bosons per lattice site as 
n,B,c = nsi^ — ns)- Here, the number of bosons per lat- 
tice site \s ub ~ n/2 — F . Therefore, in the dilute limit 
when ^ 0, almost all of the bosons are condensed 
such that ub.c ~^ ^b, however, in the dense limit when 
ub — > 1, almost all of the bosons are non-condensed 
such that nB,c ^ ^b- Making the identification that 
the number of condensed bosons is ub.c — — F), 
we obtain Ubb = 2gal as the repulsive particle (hole) 
boson-boson interaction such that the Bogoliubov rela- 
tion Vg = UBBnB.c/ {nT-Bdc) is recovered. 

Having discussed Gaussian fluctuations at zero tem- 
perature, we analyze next Gaussian fluctuations near 
the critical temperature and derive the time-dependent 
Ginzburg-Landau (TDGL) functional. 



B. Gaussian Fluctuations 
Near the Critical Temperature 

In this section, we present the results of finite temper- 
ature Gaussian fluctuations near the superfluid critical 
temperature Tc. We define the field A(x) to be the fluc- 
tuation around the order parameter saddle point value 
|Ao| = 0, and use the same procedure as in the homoge- 
neous case 0, to obtain the TDGL functional 



E 

<i,j> 



K{x) = 0, (19) 



in the position and time x = (r, t) representation. The 
expansion coefficients a, 6, q j and d are given in App. |B] 



Notice that Ci 



cSi,j is isotropic for the s-wave inter- 



actions considered in this manuscript, where Sij is the 
Kronecker delta. Next, we discuss analytically the TDGL 
functional in the weak and strong attraction limits. 



In the weak attraction (BCS) limit when < Jig- < 
AVtg, it is difficult to derive analytical expressions in 
the presence of population imbalance, and thus we limit 
our discussion only to the P = limit. In this case, 
we obtain a = Af{0) ln(T/Tc) for the coefficient of the 
linear term in A{x), b = 7Af{0)C{3)/{8Tr^T^) for the cof- 
ficient of the cubic term, c = 7J\ft{0)C{3)/{iVTT^T^) for 
the coefficent of the operator V^, and d = M{0){2'Dt+ — 
Jl+)/[2]l+{4Vt+-Jl+)]+iTrU{0)/{8Tc)] for the coefficient 
of the operator d/dt. Here, J\ft{x) = I]k J(^»^k,T)^ + 
(Vi^k,x)^]'^(^k,+ — 2^) represents an effective 'kinetic' den- 
sity of states, and the parameter Tc = (2/7r)[(2I?t+)^ — 
(m+ - 2I?i+)2]i/2exp[7 - l/(gAA(0))] is the BCS criti- 
cal temperature where 7 = 0.577 is the Euler's constant, 
= F/Af{0), and ((x) is the Zeta function. Notice that 
Tc is quite small for weak attractions, but increases with 
growing values of g. In addition, notice that the imagi- 
nary part of d is large, indicating that Cooper pairs have 
a finite lifetime and that they decay into the continuum 
of two-particle states. 

In the strong attraction (BEC) limit when w 
|efc|(l — Pe)/2 » 2Pi+, we obtain a = ai + 02 = 

-[2/1+ - Ebil - Pe)]/[gH^ - Pe)] + PelW - Pe)] for the 

coefficient of the linear term in A(a;), 6 = 61 -I- 62 = 
2/[g^{l^Pcf] - {dpc/d'ilc)/[9''{l~Pe)] for the coflicient 
of the cubic term, c — 4a^t|t|/[5'^(l — pe)^] for the co- 
efficent of the operator V^, and d = l/[g'^{l — Pe)] for 
the coefficient of the operator d/dt. Here, th = —g is the 
binding energy, e (— e) labels the excess (non-excess) type 
of fermions, and Pc — |n| — n J is the number of unpaired 
fermions per lattice site. Notice that the imaginary part 
of d vanishes in this limit, reflecting the presence of long 
lived bound states. 

Through the rescaling '^{x) = \/dA{x), we obtain the 
equation of motion for a mixture of bound pairs (molec- 
ular bosons) and unpaired (excess) fermions 



AiB«'(a;) + [UBB\-9{x)f + UBFPeix)]-i'{x) 



V2*(x) d-9{x) 

— 7 

2mB 



dt 



-0, 



(20) 



with boson chemical potential pB = —ai/d = 2/i+ — 
£(,(1 —Pe), mass mB = d/c = g/{Aa^^ti), and repulsive 
boson-boson Ubb — hia^/d^ — 2ga^c ^'^d boson- fermion 
Ubf — 0-20-c/ (dpe) — go-c interactions. Notice that, the 
repulsion Ubb between bound pairs is two times larger 
than the repulsion Ubf between a bound pair and a 
fermion, reflecting the Pauli exclusion principle. Fur- 
thermore, Ubb and Ubf are strongly repulsive due to 
the important role played by the Pauli exclusion prin- 
ciple in the lattice [33, in contrast to the homogenous 
case 0, [H, m, m where Ubb and Ub f are weakly re- 
pulsive. 

This procedure also yields Pe{x) = [a2/d + 
b2\'i'{x)\'^ /d^]/UBF, which is the spatial density of un- 
paired fermions 

Peix)^Pc-galidpc/d^ic)il-Pe)\-^ix)\^>0. (21) 
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The critical temperature in the case of zero popula- 
tion imbalance can be obtained directly from the effec- 
tive boson mass ms, using the standard BEC condition 

Tbec - 27r[nsC(3/2)]^/V"iB leading to oc t^tjg, 
which decreases with growing attraction g. Here, ^(x) 
is the zeta function. Notice that Tc in the lattice case 
vanishes for 17 ^ cx3 unlike the homogeneous case which 
saturates at Tc ~ 0.218eF, where ep is the Fermi en- 
ergy 0| • The decrease of Tc with increasing g in the 
BEC regime, combined with the increase of Tc with in- 
creasing g in the BCS regime leads to a maximum in 
between as already noted in the literature d, [sO] • 

Further insight into the differences between the homo- 
geneous and the lattice systems can be gained by com- 
paring the Ginzburg-Landau coherence length £,{T) = 
[c/(2a)]^/^, which in the vicinity of Tc becomes £,{T) = 
\gl[TcI{Tc - T)]i/2. The prefactor of the coherence 
length is 



c 
27; 



da 
df 



1/2 



(22) 



Tc. While Cgl = 

is large in the weak at- 



and is evaluated at 

[7M(0)C(3)/(8I?^2y2)]l/2 

traction limit in comparison to the lattice spacing Uc, 
S,GL — dcit-^ti/Fy^^ /g is small and phase coherence is 
lost gradually as g increases in the strong attraction 
limit when g/t+ ^ 1. This latter result is in sharp 
contrast with the homogenous case Q, where ^gl is 
large compared to interparticle spacing in both BCS 
and BEC limits, and it has a minimum near pt_|_ « [3]. 
However, ^ql depends inversely on the square-root of 
the density in both homogenous and lattice systems. 

Since the boson-boson and boson-fermion interactions 
are strongly repulsive in the BEC limit (strong attraction 
regime for fermions) due to the important role played 
by the Pauli exclusion principle in the lattice [S^, it is 
necessary to investigate this further, as discussed next. 



STRONG ATTRACTION (MOLECULAR) 
LIMIT 



A. Effective Lattice Bose-Fermi Action 

In the limit of strong attractions between fermions 
{g/t+ ^ 1), we obtain an effective Bose-Fermi lattice 
action 13511 



C-eff 



dT 



BF 



(23) 



where H'j^p = Kp 



Kb -\- Hbf 
part of the excess fermions is 



Hbb- Here, the kinetic 



(24) 



the kinetic part of the molecular bosons is 



Kb = -MB ^l^^ - *B 5Z ^\^o\ (25) 

the interaction between molecular bosons and excess 
fermions is 

Hbf = Ubf J2 flf^^l^^ + ^bf J2 flf^'^h ^ (^6) 

and the interaction between two molecular bosons is 
Hbb = Ubb J2 ^''MI^^ + Vbb J2 ^M^r (27) 

The total number of fermions is fixed by the con- 
straint n = 2nB + Pe, where ub — Nb/M is the num- 
ber of bosons per lattice site. The important parame- 
ters of this effective Hamiltonian are the excess fermion 
transfer energy tp = tc, the molecular boson transfer 
energy Ib — 2t^ti/g, the boson-fermion effective repul- 
sion Ubf — g and the boson-boson effective repulsion 
Ubb = 25. Notice that, on-site interactions Ubb and 
Ubf become infinite (hard-core) when g —t 00 as & man- 
ifestation of the Pauli exclusion principle. In addition, 
there are weak and repulsive nearest neighbor boson- 
boson Vbb oc (i| -I- i|)/5 and boson-fermion Vbf oc t^/g 
interactions. These repulsive interactions in optical lat- 
tices lead to several insulating phases, depending on 
fermion filling fractions, as discussed next. 



In the strong attraction limit, the system can be de- 
scribed by an action containing molecular bosons and 
excess fermions [l^ [2^ [2^. The existence of an opti- 
cal lattice produces different physics from the homoge- 
neous case in the strong attraction limit, because of the 
strong repulsive interaction between molecular bosons 
and excess fermions [s^. The main effect of these 
strong effective interactions is the emergence of insulating 
phases from the resulting Bose-Fermi mixture of molec- 
ular bosons and excess fermions in a lattice. 



B. Emergence of Insulating Phases 

When the fermion attraction is sufficiently strong, the 
lines BC and CD must describe insulators, as molecular 
bosons and excess fermions are strongly repulsive. In the 
following analysis, we discuss only two high symmetry 
cases; (a) nij — — I; and (b) nij = 1 and = 1/2, or 
n| — 1/2 and = 1. 

Case (a) is indicated as point C in Figs. [1] [2] and [3] 
where ni — ni = 1 such that Pe — 0. This symmetry 
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FIG. 4: (Color online) n| versus rij phase diagram in the 
strong attraction limit, indicating the CSE (superfluid), the 
metallic (lines AB and DE), and insulating (lines BC and CD) 
phases with special points Fermi-Pauli insulator (FPI) , Bose- 
Fermi checkerboard (BFC), Bose-Mott insulator (BMI), and 
BMI/FPI phase separation. 



point is a Fermi-Pauli (band) insulator for weak attrac- 
tion since both fermion bands are fully occupied, and 
a Bose-Mott Insulator (BMI) in the strong attraction 
limit. There is exactly one molecular boson (consisting 
of a pair of t and | fermions) at each lattice site, which 
has a strong repulsive on-site interaction with any addi- 
tional molecular boson due to the Pauli exclusion prin- 
ciple. Notice that this case is very similar to the atomic 
BMI transition observed with one-species atomic Bose 
systems 29]. 

In case (a) , H'^p reduces to a molecular Bose-Hubbard 
Hamiltonian with the molecular Bose filling fraction 
IT'S = n'/2 — F, thus leading to a molecular BMI when 
ub — 1 beyond a critical value of Ubb- A schematic 
diagram of this phase is shown in Fig. EJ^b). The crit- 
ical value Ugg needed to attain the BMI phase can 
be estimated using the approach of Ref. [i^l leading to 
U 



3(3-1- V8)tB, which in terms of the underlying 
fermion parameters leads to gc = A.lSy/titi for the crit- 
ical fermion interaction. This value of gc is just a lower 
bound of the superfluid-to-insulator (SI) transition, since 
H^^p is only valid in the g ^ t+ limit. Some signatures 
of this SI transition at g,, have been observed in recent 
experiments p^ . 
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FIG. 5; (Color online) Schematic diagrams for the (a) Fermi- 
Pauli insulator (FPI), and (b) Bose-Mott insulator (BMI) 
phases. 

We would like to make some remarks on the validity 



Hgp presented above, which intrinsically assumes that 
the molecular bosons arc small in comparison to the lat- 
tice spacing. A measure of the 'smallness' of these molec- 
ular bosons is the average size of the fermion pairs defined 
(for pe = 0) by 



,pair 



(0(k)|r2|0(k)) 



1/2 



(28) 



where (/)(k) = Ak/(2i?k,+) is the ground state pair wave 
function. While S,pair = vf/{^\^o\) is large in the weak 
attraction limit in comparison to the lattice spacing Oc, 
Cpair = 2\/2Pac|l — 2F\t+/g in the strong attraction 
limit when g/t^ ^ 1. These expressions are valid for 
low and high filling fractions, F ^ and ^ 1, re- 
spectively. At g — gc, we obtain ^. 



pair 



l.lTflc for 

i| = and ^pair/o.c ~ 1.74ac for = O.lSij. Notice 
that larger values of g lead to molecular boson sizes that 
are smaller than the lattice spacing, thus validating the 
effective Bose-Fermi description derived above for g ^ t+ 
to values of g close to gc- In the lattice case, the results 
above suggest that even if the size of the bound pairs 
decreases as t^ad g with increasing g, the Pauli pressure 
prohibits two bound pairs or a bound pair and an excess 
fermion to occupy the same lattice site, in the single- 
band description discussed here. This indicates that the 
bound pairs do not loose their fermionic nature and can 
not be thought as structureless in lattices, which should 
be contrasted with the homogenous case where bound 
pairs become structureless with increasing attraction. 

Case (b) is indicated as crosses in Figs. [U [2] and |4] at 
points n| = l,nx = 1/2 or — 1/2, ni — 1, where the 
molecular boson filling fraction ns = 1/2 and the ex- 
cess fermion filling fraction is Pe = 1/2. At these high 
symmetry points, molecular bosons and excess fermions 
tend to segregate, either producing a domain wall type 
of phase separation between a molecular Bose-Mott in- 
sulator (BMI) region and a Fermi-Pauli insulator (FPI) 
region or a checkerboard phase of alternating molecular 
bosons and excess fermions (BFC) depending on the ra- 
tio Vbb/Vbf- a schematic diagram of these two phases 
is shown in Fig. [6] 
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FIG. 6; (Color online) Schematic diagrams for the (a) Bose- 
Fermi checkerboard (BFC) phase, and (b) Bose-Mott insula- 
tor (BMI) and Fermi-Pauli insulator (FPI) phase separation. 

The checkerboard phase shown in Fig [6l^a) is favored 
when Vbb > 2Vbf- At the current level of approxi- 
mation, we find that when t-^ — ti phase separation is 
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always favored, however when | (J,) fermions are in ex- 
cess the checkerboard phase is favored when > VS^t 
{ti < <|/\/3). Therefore, phase separation and checker- 
board phases are achievable if the tunneling ratio rj can be 
controlled experimentally in optical lattices. Notice that 
this checkerboard phase present in the lattice case is com- 
pletely absent in homogeneous or harmonically trapped 
systems [H [H [H i2S] . 

Thus, the strong attraction limit in optical lattices 
brings additional physics not captured at the saddle 
point, and not present in homogenous or purely har- 
monically trapped systems. Having discussed the strong 
attraction (molecular) limit and its possible phases, we 
comment next on their experimental detection. 



VI. DETECTION OF INSULATING PHASES 

It may be very difficult to detect all the superfluid and 
insulating phases proposed in situ, and thus one may 
have to turn off the optical lattice to perform time-of- 
flight measurements. The phases proposed here in the 
molecular limit could be detected either via the measure- 
ment of momentum distribution, density-density corre- 
lations or density Fourier transform. For instance, the 
Fermi-Pauli insulator (FPI), Bose-Mott insulator (BMI), 
Bose-Fermi checkerboard (BFC), and phase separated 
BMI/FPI do not exhibit phase coherence in their mo- 
mentum distribution, while the superfluid phase does. 
Therefore, through the momentum distribution measure- 
ments, one should be able to differentiate between super- 
fluid and insulating phases, which has been possible for 
Bose systems in optical lattices [2^. However, in order 
to detect clearly the different insulating phases, it may 
also be necessary to explore the Fourier transforms of the 
densities, or of the density-density correlation functions. 

The relevant density-density correlation functions to 
characterize the insulating phases are CF,F(r, r') — 
(rii?(r)rii?(r')) which correlates the excess fermions, 
C_B,s(r, r') — (nB(r)ns(r')) which correlates the molecu- 
lar bosons, and CB.F(r,r') = ("-s(r)nF(r')) which corre- 
lates the molecular bosons and excess fermions. Defining 
the relative R = r — r' and the sum Rc = (r + r')/2 of 
the coordinates, we can write the correlation functions 
as Ca,/3(R, Rc) = {riaCRc + R/2)n/3(Rc - R/2)), where 
a and /3 are {F, B}. This leads to the Fourier transform 
Ca,/3(q, Qc), where q is the relative and Qc is the sum of 
the momentums. For a translationally invariant system 
considered here, it is sufficient to look only to qc = 0, and 
define Ca,i3{q) = ^^^^(q, 0). Therefore, for the FPI, BMI 
and BFC insulating phases, correlation functions have 
the generic form Ca,f3{q) = Y.i,j exp[-iq • (n^^ - rj^^)], 
where and rj ^ are the locations of a and f3 particles. 

For instance, in the case of a two-dimensional system, 
the FPI phase [shown in Fig. [DJa)] has a strong peak in 
the correlation function Ci?_i?(q) at q = (0,0). Similarly, 
the BMI phase [shown in Fig. [5t^b)] has a strong peak in 



the correlation function Cb.b{(1) at q = (0, 0). However, 
the BFC phase [shown in Fig. [IKa)] has strong peaks in 
the correlation functions CF._F(q) and CB.s(q) both at 
q — (0,0) and q = ±(7r,±7r), while it has strong peaks 
in the correlation function CB,F(q) at q = (0,0). 

The insulating FPI, BMI, and BFC phases are modi- 
fied in the presence of an overall harmonic trapping po- 
tential, since the molecular bosons experience the har- 
monic potential Vsir) — {m^uj'^ + TOj^(jjJ)r^/2, and the 
excess fermions of type a experience a different har- 
monic potential Vp.air) — maLo^r^ jl. Provided that the 
harmonic potentials are not too strong, the character- 
istic peak locations of the correlation functions remain 
essentially the same, but additional broadening due to 
the inhomogeneous nature of the trapping potential oc- 
curs. However, when the harmonic potentials are sufR- 
ciently large, important modifications of these phases are 
present, including the emergence of complex shell struc- 
tures. 

Lastly, the BMI/FPI phase separated state is strongly 
modified even in the presence of weak trapping poten- 
tials. The single domain wall density structure [shown in 
Fig. [Sfb)] turns into a new density profile with a region 
where the bosons tend to stay at the center of the trap- 
ping potential due in part to their tighter confinement, 
and the excess fermions are pushed out of the center by 
the strong Bose-Fermi repulsion. The fermions then sur- 
round the bosons, and standard in situ radio-frequency 
measurements can detect the existence of the BMI/FPI 
phase separated state, as done recently with population 
imbalanced fermion mixtures [l|, 0, H, 3 ■ 

Having commented briefly on the experimental detec- 
tion of the superfluid and insulating phases, we present 
our conclusions next. 



VII. CONCLUSIONS 

Using an attractive Fermi-Hubbard Hamiltonian to de- 
scribe mixtures of one- or two-species of atoms in op- 
tical lattices, we obtained the ground state phase dia- 
gram of Fermi- Fermi mixtures containing normal, phase- 
separated and coexisting superfluid/excess fermions, and 
insulating regions. We discussed the cases of balanced 
and imbalanced populations of Fermi-Fermi mixtures 
such as ^Li or only; and mixtures of ^Li and ^°K; ^Li 
and *^Sr; or "^^K and *^Sr. We showed that population 
imbalanced Fermi-Fermi mixtures reduce to strongly in- 
teracting (repulsive) Bose-Fermi mixtures in the molecu- 
lar limit, in sharp contrast to homogenous systems where 
the resulting Bose-Fermi mixtures are weakly interacting. 
This result is a direct manifestation of the Pauli exclusion 
principle in the lattice case, since each Bose molecule con- 
sists of two fermions, and more than one identical fermion 
on the same lattice site is not allowed, within a single- 
band description. This effect together with the Hartree 
energy shift lead to a flUing dependent condensate frac- 
tion and to sound velocities which do not approach zero, 
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in contrast to the homogenous case. 

Furthermore, we showed that several insulating phases 
appear in the strong attraction limit depending on fill- 
ing fraction and population imbalance. For instance, we 
found a molecular Bose-Mott insulator (superfluid) when 
the molecular filling fraction is equal to (less than) one for 
a population balanced system where the fermion filling 
fractions are identical. When the filling fraction of one 
type of fermion is one and the filling fraction of the other 
is one-half (corresponding to molecular boson and excess 
fermion filling fractions of one- half), we also found ei- 
ther a phase-separated state consisting of a Fermi-Pauli 
insulator (FPI) of the excess fermions and a molecular 
Bose-Mott insulator (BMI) or a Bose- Fermi checkerboard 
(BFC) phase depending on the tunneling anisotropy ra- 
tio. 

All of these additional phases and the possibility of ob- 
serving more exotic superfluid phases in optical lattices 
with p-wave order parameters [3l| make the physics of 
Fermi-Fermi mixtures much richer than those of atomic 
bosons or Bose-Fermi mixtures in optical lattices, and 
of harmonically trapped fermions. Lastly, the molecu- 
lar BMI phase discussed here has been preliminarily ob- 
served in a very recent experiment J.2,], opening up the 
experimental exploration of the rich phase diagram of 
fermion mixtures in optical lattices in the near future. 

We thank National Science Foundation (DMR- 
0709584) for support. 



corresponding to the |qp term, and 



(A3) 



corresponding to the w'^ term. 

The coefficients necessary to obtain the diagonal 
phase-phase matrix element are 

Q - §E{^'^°^(^-«^)-[«+-3|AoP)t^ 



corresponding to the |qp term, and 



(A4) 



(A5) 



corresponding to the vo^ term. 

The coefficient necessary to obtain the off-diagonal ma- 
trix element is 



1 V- CkH 



(A6) 



corresponding to the w term. These coefficients can be 
evaluated analytically in the BCS and BEC limits, and 
are given in Sec. IIV Al 



APPENDIX A: EXPANSION COEFFICIENTS AT 
ZERO TEMPERATURE 



APPENDIX B: EXPANSION COEFFICIENTS 
NEAR THE CRITICAL TEMPERATURE 



In this Appendix, we perform a small q and w ex- 
pansion of the effective action at zero temperature (T = 
0) [1^. In the amplitude-phase basis, we obtain the ex- 
pansion coefficients necessary to calculate the collective 
modes, as discussed in Section IIV Al We calculate the 
coefficients only for the case of s-wave pairing with zero 
population imbalance P = 0, as extra care is needed 
when P 7^ due to Landau damping. In the long- 
wavelength (|q| — > 0), and low- frequency (w — s- 0) limits 
the condition {ii;,eq,+} ^ min{2£'k,+} is used. 

The coefficients necessary to obtain the diagonal 
amplitude-amplitude matrix element are 



A- 



1 



|Ao| 



M ^ 2El^ 



(Al) 



corresponding to the (q = 0, w = 0) term, 
, f £;^,+ -3 |Aop 



^ = M^X 4^5 cos(/c,ac) 



10|Ao|^ 



(A2) 



In this Appendix, we derive the coefficients a, 6, Ci_j 
and d of the time dependent Ginzburg-Landau theory 
described in Section IIVBI Wc perform a small q and 
w expansion of the effective action near the critical tem- 
perature (T « Tc), where we assume that the fluctuation 
field A(r, t) is a slowly varying function of position r and 
time t |28l|. 

The zeroth order coefficient is given by 



1 1 ^ 

a = > 



|rs(k)|^ 



(Bl) 



where X\s_^± = (-'^k.T i ^Kl)/'^ ^'^d Xk.^ 
tanh(/3^k,a/2). 

The second order coefficient is given by 



M 



E { (^k,T>k,T<T + ^k,i>k,itx) sin(/c,ac; 

k 

0^ [4t_ sin(fciac) sin(fcjac)C_ 



o4k, 

— 2(5i J- cos(fcjac)C+ 



At^_ sin(fciac) sin(A;jac) 



k,-l- 



?kH 



-'^k,+ 



}|rs(k)p, (B2) 
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where C± = (Fk.T^T ± ^k4ii)/2, and Fk.o 
sech^(/3^k,(T/2). Here, Sij is the Kronecker delta. 
The fourth order coefRcient is given by 



-E 



|r.(k)r 



(B3) 



The time-dependent coeificient has real and imaginary 



parts, and for the s-wave case is given by 



= lim y 



M 



4e 



|r.(k)(i34) 



where S{x) is the delta function. These coefficients can 
be evaluated analytically in the BCS and BEC limits, 
and are given in Sec. IIVBI 
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